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Abstract 

We study the Bose-Einstein correlations of two photons and their coherent 
properties that can provide the information about the space-time structure 
of the emitting source through the Higgs-boson decays into two photons. We 
argue that such an investigation could probe the Higgs-boson mass. The 
model is rather sensitive to the temperature of the environment and to the 
external distortion effect in medium. 

PACS numbers: 14.80 Bn, ll.lO.Wx 



Introduction.- The hadron colliders Tevatron and the LHC are at the stage to 
provide particle physicists with treasure of data. These data allowed a precise mea- 
surement of many important parameters of particle physics in order to test their 
consistency, and to discover new particles, in particular, the Higgs-boson. Today 
little is known about the interplay between the Higgs and gauge bosons in medium. 
The method of Bosc-Einstcin (BE) correlations is clear and undeniable part of this 
physics, complicating the quantum statistical description of multiparticle final states 
(for review of BE correlations see [1]). 

Historically, the BE correlations measurements have concentrated on pion pairs 
correlations, however have also been applied to more heavy hadrons and gauge 
bosons. 

Correlations of two photons can provide information about the space-time struc- 
ture of hot matter prior to freeze-out. However, the BE correlations with direct 
(primary) photons is faced to difficulties compared to hadrons primarily due to the 
small yield of photons emitted directly from the hot region just after particle colli- 
sions. The main background is related with decays of hadrons into photons. In order 
to reject most of the hadron background, the reconstructed events in a spectrometer 
were required to have a deposit energy (level) of greater than the energy well above 



the minimal energy at which photons emerge from the hadrons. Experimentally, the 
background constitutes a major difficulty which, however, to some extent, can 
be taken care of by measuring the photon pair invariant mass. 

In Higgs searching experimental program, it is supposed that two photons are 
produced mainly through the decay of light Higgs-boson H. The H ^ '-ff decay 
mode is one of the promising discovery channels for the Standard Model (SM) Higgs- 
boson in low mass region (114 < Mh < 150 GeV). Despite the small branching ratio 
(~ 10^^ for Mh = 120 GeV) this channel has a simple signature and a very good 
mass resolution 1.5 GeV). From the theoretical point of view the decay of the 
Higgs into two photons emerges through the one-loop composed of either fermions 
(e.g., quarks) or vector bosons, or the scalar particles (see, e.g., [2] and the references 
therein) . 

The two-photon BE correlations in Higgs physics can provide a powerful tool to 
estimate the Higgs-boson mass via the characteristic stochastic scale Lgt, related to 
correlation radius, which defines the geometrical size of the two-photon source. The 
relative correlation function is strongly dependent on the space-time properties of 
the Higgs decays. 

We suppose the gamma-quanta are created by some random currents (sources) 
which exist in a restricted space-time region and they are chaotically and randomly 
disturbed by external fields (forces). The restricted source area is characterized by 
the scale Lgt, the meaning of which is explained in [3]. 

In the quark-loop scheme, the non-relativistic bremsstrahlung formula for the 
current of one photon emission (with four-momentum p — {p^, p)) is 

mqp^ 

where q is the difference between the spatial momenta of a quark and an antiquark 
with the mass m^, e\{p) is the vector of linear polarization of a photon; g is the 
coupling constant; Eq means the phenomenological parameter which depends on the 
initial energy. For two sources the transition current is 

2 

j\p) = Y.^'''-3n{p)- 

n=l 

Index n labels the independent quark-antiquark interactions taking place at different 
space-time points yn which are considered to be randomly distributed in the region 

of the (photon) source. 

By studying BE correlations of identical particles, it is possible to determine the 
time scale and spatial region over which particles do not have the interactions. Such 
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a surface is called as decoupling one. In fact, for an evolving system, such as one 
for pp collisions, it is not really a surface, since at each time there is a spread out 
surface due to fluctuations in the last interactions, and the shape of this surface 
evolves even in time. 

Let us consider the distribution function 

Nl2iPl,P2) = {\j''iPl)j'^iP2)\') (1) 

normalized to the product iVi(pi) • A^2(p2) with one-particle distribution functions 
Ni{pi) = {\J^'{Pi)\'^) {i = 1,2). The function ([l| defines the probability to find 
two photons with momenta pi and p2 issued at the points yi and y2- The crossing 
momenta has to be taken into account. 

In this paper we suggest the proposal for the measurements of quantum corre- 
lations between two photons at hadron colliders. This is theoretically supported by 
the quantum field theory model approaches at finite temperature [4,3], where one 
of the main parameters is the temperature of the particle emitting source. 

We imply that photons do not strongly interact with the medium: they carry 
information about early stage of reaction. Our pragmatic definition is: photons are 
produced not from hadronic decays. The gamma-quantum can be measured from the 
virtual yield which is observed as some mass of the spectrum for lepton-antilepton 
pair. The main channels are: pp ^ H ^ 7*7* — )■ 2/i~2/i"'", 2e~2e~^, e~ fi~ , ... 
in, e.g., pp collisions. 

BE correlation function. - In general, a pair of identical bosons with momenta 
Pi and P2 and the mass m produced incoherently from an extended source will have 
an enhanced probability 6*2(^1,^2) = Ni2{pi, P2) / [Ni{pi) ■ A^2(P2)] to be measured in 
terms of differential cross section cr, where 

AT / \ 1 ^^"^ 
iVl2 Pi,P2 = - 

a a\li d\l2 

to be found close in 4-momentum space 3^4 when detected simultaneously, as com- 
pared to if they are detected separately with 

1 da _ d?pi 

ad^i' ' ~ (277)3 2Ep, 



MP^) = -^^ = ^ , Ep^ = J0f + m\ 1 = 1, 2. 



In an experiment, one can account the inclusive density ^2(^1,^2) which describes 
the distribution of two particles in f2 (the sub- volume of the phase space) irrespective 
of the presence of any other particles 

11 d'^n2 

P2{PUP2) 



2! Uevents dpi dp2 
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where n2 is the number of particles counted in a phase space region {pi+dpi, p2+dp2) ■ 
The multiphcity normahzations stand as 

p{p)dp={N), 



n 



P2{pi,P2) dpi dp2 = {N{N - 612)), 



where (A^) is the averaged number of produced particles (mean multiplicity). Here, 
S12 = for different particles, while 5i2 = 1 in case of identical particles (coming 
from the same event). 

The following relation can be used to retrieve the BE correlation function C2{Q)'- 

C (O) - ^^^^ (2) 

where N{Q), in general case, is the number of particle pairs (off-shell photons) in 
BE correlation pattern with 

Q = ^J-{pi-P2)^-{pi-P2Y = VM^-Am^. (3) 

In definitions ([2| and ([3|, A ^'''^-^ is the number of particle pairs without BE correla- 
tions and M = \Jipi + ^2)^ is the invariant mass of the pair of photons. 

An essential problem is the estimation of the reference distribution N"^^^ {Q) in 
([2]). If there are other correlations besides the BE ones, N^'^^{Q) should be replaced 
by a reference sample corresponding to two-particle distribution in a geometry with- 
out BE correlations. Hence, the expression ^ represents the ratio between the 
number of 7*7* pairs N{Q) in the real world and the reference sample N'^^^ {Q) in 
the imaginary world. Note, that N'''^^ can not be directly observed in an experiment. 
Different methods are usually applied for the construction of the reference samples 
[1] , however all of them have strong restrictions. One of the preferable methods is to 
construct the reference samples directly from data. The 7*7* correlations can be es- 
timated for each bin of the photon average transverse momentum pt = \pti + PT2 1 / 2 
as the ratio of the distribution of photon pair invariant relative momenta where both 
photons with transverse momenta and Pt2 were taken from the same event to 
the same distribution but with the photons of the pairs taken from different (mixed) 
events. 

In general, the shape of C2{Q) is model dependent. The most simple form for 
C2{Q) is often used for experimental data fitting [5] 

C2iQ) = Co-[l + Ae-(«-^)1 ■ (1 + eQ), (4) 
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where Cq is the normahzation factor, A is the coherent strength factor, meaning 
A = 1 for fully incoherent and A = for fully coherent sources; the symbol R is 
often called as the " correlation radius" , associated with the measure of the particle 
source size. The linear term in Q is supposed to be account within the long-range 
correlations outside the region of BEC. For the experimental data fitting it is often 
used either k = 1 or k, = 2. It is assumed that the maximum of C2{Q) is 2 for pi = p2 
if no any distortion and final state interactions are taking into account. Note that 
the distribution of bosons can be either far from isotropic, usually concentrated in 
some directions, or almost isotropic, and what is important that in both cases the 
particles are under the random chaotic interactions caused by other fields in thermal 
medium. In the parameterization Q all of these issues should be embedded in A. 

However, it is known that measured both the size of the source and A are de- 
pendent on the transverse momentum pt (or the transverse mass m^) and the 
multiplicity [6]. The model which predicts the pt (mj-)-, N- and the temperature 
T - dependencies as well as the stochastic forces influence is the so-called Stochas- 
tic Evolution Model at finite temperature, SEM13, [3]. In this article we develop 
this model further and apply it to reconstruction of the two-photon correlations. 
The two-particle correlation function in SEM^ is (the correlations at large Q are 
neglected) 

C2{PuP2) ^ i{N) {1 + Ai(/3) e-^-« [1 + A2(/3) e^-^/^j | ^ 

where Aps^ = —{pi —P2Y ^f^v (Pi ~P2Y is the smearing smooth generalized function. 
3fJ^j, is the nonlocal structure tensor of the space-time size and it defines the region 
of emitted photons. ^(A^) depends on as i{N) = {N{N - 1))/{N)'^. The function 
Ai(/3) is the measure of the strength of BE correlation between two photons, Ai(/3) ~ 
7(0;, /?)/(! + a)^, and the correction is A2(/3) ~ 2al ^^{uo, /3). The function 7(0;, /3) 
calls for quantum thermal features of BE correlation pattern and is defined as 

7(u;,/3)^7(n) = ^^^^ n(c.O ' "^^^ " = ^ = " ' 

where n{u,(3) is the mean value of quantum numbers for Bose-Einstein statistics 
particles with the energy u in the thermal bath with statistical equilibrium at the 
temperature T = 1/13. The following condition 'Y^-nj{uj, fi) = N is evident, where 
the discrete index j stands for the one-particle state. 

The function a = a{/3,m) entering Ai and A2 in ([s]), the measure of stochastic 
forces, summarizes our knowledge of other than space-time characteristics of the 
particle emitting source, and it varies from to 00. 
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In terms of time-like Lq, longitudinal Li and transverse Lt components of the 
space-time size -^/T^, the distribution ApSR in ^ looks like 



pL 



(7) 



where Lq is treated as the measure of the particle emission time, or even it represents 
the interaction strength of outgoing particles. 

Hence, we have introduced a new parameter y^^^, which defines the region of 
nonvanishing particle density with the space-time extension of the particle emission 
source. Formula ^ must be understood in the sense that exp(— Aps^) is a distribu- 
tion that in the limit L — oo strictly becomes a 5 - function. For practical using 
with ignoring the energy-momentum dependence of a, one has 



c,{Q) ^ e(iv) 



7(n) 



+ 



a] 



where the stochastic scale Lgt = Lst{P) is the measure of the space-time overlap 
between two photons, and the physical meaning of Lgt depends on the fitting of 
C2(Q)-function. Lgt can be defined through the evaluation of the root-mean-squared 
momentum Qrms as: 
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'C2{Q,(3)-l 





where Lgt and Qrms{(3) are related to each other by means of 



1 + 



1 



1 



4a(/3) V 2 



1/2 



The following restriction a/3/2 < {Lst ■ QrmsiP)) < is evident, where the lower 
bound satisfies to the case a — ?■ (no any distortion in the particle source), while 
the upper limit is given by the very strong influence of chaotic forces, a — )■ oo. The 
result is rather stable in the wide range of variation of a. 

In-medium distortion.- It has been emphasized [3] that there are two different 
scale parameters when the BE correlations are studied. One of them is the corre- 
lation radius R introduced in (|4]). In fact, the latter gives the "pure" size of the 
particle emission source. The other parameter is the scale Lgt of the production par- 
ticle region where the stochastic, chaotic distortion due to environment (medium) 
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is enforced. Lgt carries the dependence of the particle mass, the a-coherence degree 
and the temperature. 

The question arises: how can BE correlations be used to determine the effective 
stochastic scale Lgt and, perhaps, the phase transition? We suppose the changes of 
7*7* production region and effects yielding the dynamical variables and parameters 
of BE correlations due to in-medium distortion. 

In the Higgs-boson rest frame, there is a kinematical configuration for two pairs 
of final lepton, antilepton momenta {pi,Pj) and {qi, qj): pi ~ p'j and qi ~ q'j. The final 
state mimics a two-body final state, and if the leptons being the electrons or muons, 
the virtual photons with momenta p and p' are nearly on mass-shell, p"^ c:^ p' ~ 4 mf, 
where mi is the lepton mass. For the configuration given above, the stochastic scale 
Lgt has different behavior depending on T. At lower temperatures we have 



48a(A^)mJ^/^T3/2 (l 



15 T_ 

16 mi 



(9) 



where the condition 2 n /3 > 1 is taken into account for any integer n; is the 
mass of the Higgs-boson. 

On the other hand, at higher temperatures, one has 



48C(3)a(Ar) T3 mf 



C(3) = ^n-3 = 1.202. 



(10) 



n=l 



In case of two real photons one can use the longitudinal stochastic scale L^°"^(m; = 

mT/2), with the average transverse mass m^ = 0.5(y^p|r-|- \Jp\^) in the frame of, 

e.g., the Longitudinal Center of Mass System (LCMS) [7]. 

It turns out that the scale L^t defines the range of stochastic forces. This effect 
is given by a(A^)-coherence degree which can be estimated from the experiment 
through the function C2{Q) as Q close to zero, 6*2(0) 



a{N) 



1 + ^'l\n) - 62(0) + ^"\n)^J C2mi^'^{n) -2] + 2 

6*2(0) - 1 ' 



where 7(n) is defined in ([6j) and 6*2(0) = C*2(0)/^(A^). The upper limit of 6*2(0) 
depends on (A^) and the quantum thermal features of BE correlation pattern given 
by 7(n): 6*2(0) < ^(iV)[l — 7^/^(n)/2]~^. This upper limit is restricted by the 
maximal value of 2 in the ideal case as (iV) — )• 00 and 7(n) = 1. 
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Note, that the hmit a — t- oo yields for fully coherent sources with large (A^), 
while the a case stands for fully chaotic (incoherent) sources for small (A^). 
The increasing of T leads to squeezing of the region of stochastic forces influence, 
and LstiT = Tq) = R a.t some effective temperature Tq. The higher temperatures, 
T > Tq, satisfy to more squeezing effect and at the critical temperature Tc the scale 
LstiT = Tc) takes its minimal value. Obviously, Tc ~ 0(200 GeV) defines the phase 
transition where the electroweak symmetry restoration will occur. Since in the phase 
where T > T^ all the masses tend to zero and a — )■ 0, one should expect the sharp 
expansion of the region with Lst{T > Tc) — )■ oo. 

Using the relation between L; and rriT obtained from the Heisenberg uncertainty 
relation Li(mT) = c\/h LQ/2'KmT [8] at T = Tq, one can estimate Rq = Lq within 
the formula: 

2 TT 772,7- 



Rq 



C^h 



3C(3)a(Ar)T3m^ 



(12) 



where c and h are the speed of light in vacuum and the Planck constant, respectively. 
The expression (12) relates the two-photon emission time Rq with the Higgs-boson 
mass Mh, the transverse mass of two photons mx and the temperature of two 
photons emission source Tq. We obtain Rq ~ 10~^^sec for Mh = 120 GeV with 
rriT = 1 GeV (two muon pairs in the final state) in the wide interval Tq = 10 — 100 
GeV. In Fig.l we display the Higgs-boson masses for various values of Tq and a 
with R = 1 fm and mi = = 105 MeV. As an illustration we present here an 
estimation of Mh = 138 GeV assuming Tq = 20 GeV, a = 1 (A ~ 0.8). 

The qualitative relation between R and Lgt above mentioned is the only one we 
can emphasize in order to explain the mass and the temperature dependencies of 
the source size. The dependence of the stochastic measure of chaoticity a on the 
minimum scale cut Lgt and Qrms can be used to define the fit region for different 
Pt- Such a minimum cut on Lgt introduces a lower cutoff on R. 

There are a number of effects which may give rise to Q — correlations and, thus, 
mimic the two-photons BEG coming from Higgs decays. These include i) correlations 
of light hadrons missidentified as photons, ii) radiative decays of resonances in both 
pseudoscalar and vector sectors. Hi) collective flow, etc. 

Apparatus or analysis effects which may result in the sensitivity of external 
random forces influence, may be investigated by studying the dependence of the 
correlation functions on the stochastic scale Lgf This effect is expected to contribute 
strongly at small Lgt (or large a and T). 

We find that the stochastic scale Lgt decreases with increasing T slowly at low 
temperatures, and it decreases rather abruptly when the critical temperature is 
approached. Actually, the experimental measuring of R can provide the precise 



8 




Fig. 1: Higgs-hoson masses as a function of Tq and a{N) for R = Ifm. We have 
set mi = rrin = 105 MeV. 



estimation of Mh at an effective temperature Tq. The parameter a can be extracted 
from tlie experiment using C2{Q) with Q being close to zero (11). 

Our model is consistent with the idea of the unification of weak and electro- 
magnetic interactions at T > T^, predicted by Kirzhnitz and Linde in 1972 [9]. In 
addition, there is the analogy with the asymptotic free theory: the properties of en- 
vironment (medium) are the same as those composed of free particles in the infinite 
volume (Universe). 

To summarize: we faced to the two-photon correlation function C2{Q,P) in 
which the contribution of N, T, a are presented. This differs from the methods 
based on, e.g., the Goldhaber-like approach (|4]), the trigonometrical parametrization 
[10], etc. In particular, by studying the BE correlations of two photons at finite 
temperature we can predict the mass of the SM Higgs-boson. There is the possibility 
to determine the precision with which the source size parameter and the strength 
chaoticity parameter(s) can be measured at hadron colliders. Such investigations 
provide also an opportunity for probing the temperature of the particle production 
source and the details of in-medium distortion. 

I thank I. Gorbunov for help in numerical calculations. 
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